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The  overall  objective  of  this  research  investigation  is  to 
develop  a  methodology  for  displaying  microwave  dosimetric 
data.  Special  emphasis  was  placed  on  developing  a  method 
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an  image.  A  cubic  spline  interpolation  program  was  developed 
to  uniformly  expand  the  data  set  for  presentation  on  a  large 
format  display.  A  two  segment  enhancement  curve  was  developed 
which  gives  the  operator  the  capability  of  creating  his  own 
lookup  table  for  pseudocolor  translation  of  images  stored  in 
refresh  memory.  The  results  to  date  show  that  pseudocolor  imaging 
is  an  effective  method  of  analyzing  microwave  dosimetric  data. 
Additional  processing  algorithms  need  to  be  developed  and  the 
system  hardware  should  be  optimized  for  image  processing. 
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SUMMARY 


The  overall  objective  of  this  research  investigation  is  to 
develop  a  methodology  for  displaying  microwave  dosimetric  data. 
Special  emphasis  was  placed  on  developing  a  method  for  evaluating 
the  spatial  variation  of  absorbed  microwave  energy  by  specific 
organs  and  organ  subdivisions.  Techniques  were  developed  by  a 
computer  and  displayed  on  an  image  in  pseudocolor.  Manipulation 
of  the  data  using  the  image  display  system  allows  the  operator 
to  highlight  organ  subdivisions  within  an  image.  A  cubic  spline 
interpolation  program  was  developed  to  uniformly  expand  the  data 
set  for  presentation  on  a  large  format  display.  A  two  segment 
enhancement  curve  was  developed  which  gives  the  operator  the 
capability  of  creating  his  ow,.  lookup  table  for  pseudocolor  trans 
lation  of  images  stored  in  refresh  memory.  The  results  to  date 
show  that  pseudocolor  imaging  is  an  effective  method  of  analyzing 
microwave  dosimetric  data.  Additional  processing  algorithms  need 
to  be  developed  and  the  system  hardware  should  be  optimized  for 
image  processing. 


SECTION  I 


Introduction 


The  research  described  in  this  Annual  Report  is  the  result 
of  efforts  performed  during  the  first  year  of  this  contract.  As 
a  result  of  an  inordinate  delay  in  obtaining  necessary  hardware  for 
this  project,  it  was  necessary  to  extend  the  time  for  performance 
of  the  contract  beyond  the  original  one  year.  The  overall  objec¬ 
tives  of  this  research  investigation  is  to  develop  a  methodology 
for  displaying  microwave  dosimetric  data  on  a  CRT  type  display 
for  use  in  data  interpretation  and  analysis.  A  color  image  pro¬ 
cessing  system  was  procured  and  interfaced  to  the  laboratory 
Hewlett-Packard  computer  system.  Data  obtained  from  the  microwave 
scanner  is  processed  by  the  computer  and  then  displayed  on  the 
image  system.  Using  the  image  system,  the  experimenter  can  mani¬ 
pulate  an  image  to  enhance  areas  of  interest  for  better  analysis. 
Shortly  after  this  system  was  placed  in  operation,  it  became 
apparent  that  viewing  the  data  set  in  a  one-to-one  aspect  ratio 
produced  an  image  which  was  too  small  for  detailed  analysis.  A 
cubic  spline  interpolation  program  was  developed  to  uniformly 
expand  the  image  to  a  larger  format.  Also,  a  two  segment  enhance¬ 
ment  curve  was  developed  to  allow  the  operator  to  define  his  own 
video  lookup  table  for  pseudocolor  translation  of  images  stored  in 
refresh  memory. 


SECTION  II 


DEVELOPMENT  OF  DATA  ANALYSIS  SYSTEM 


This  section  of  the  report  provides  an  overview  of  the  dis¬ 
play  system  and  a  description  of  its  application  to  data  analysis 
in  pictorial  format.  The  operation  of  the  image  display  system  is 
straight  forward  and  does  not  involve  any  computer  programming  on 
the  part  of  the  user.  A  brief  description  of  the  commands  used  by 
the  operator  to  interact  with  the  system  is  presented. 


A.  Selection  of  a  Display  Media 

The  objective  of  this  phase  was  to  select  a  method  for  pre¬ 
senting  experimental  results  in  a  manner  that  maximizes  the  infor¬ 
mation  content  and  interpretability  of  the  data.  To  date,  the 
data  display  has  primarily  involved  an  isometric  presentation  that 
represents  energy  absorption  or  phase  shift  as  a  height  above  a 
planar  surface  whose  perspective  can  be  varied.  The  major  disad¬ 
vantage  of  this  approach  is  that  the  aerial  perspective  contains 
no  information.  It  was  felt  that  an  aerial  perspective  was  criti¬ 
cal  to  the  successful  analysis  of  the  data  because  an  aerial  per¬ 
spective  would  allow  a  determination  of  the  spatial  relation  be¬ 
tween  the  organ  subdivisions. 

The  data  obtained  from  the  microwave  scanner  system  is  in  a 
digital  format.  Each  data  point  represents  the  energy  absorption 
of  a  specific  physical  location  on  the  organ.  Using  the  co-ordi¬ 
nate  information  associated  with  each  intensity  value,  the  data 
can  be  assembled  to  form  a  mosiac  "picture"  of  the  energy  absorp¬ 
tion  of  the  organ.  This  process  is  identical  to  that  employed  in 
reconstructing  an  image  with  data  from  a  satellite.  Hence,  a 
literature  search  of  commercial  image  processing  systems  was  under¬ 
taken.  Systems  from  the  following  venders  were  considered: 

.  Ramtek  Corporation,  Sunnyvale  CA 
.  Lexidata  Corporation,  Burlington  MA 
Grinnel  Systems,  San  Jose  CA 
.  Comtal  Corporation,  Altadena  CA 

The  Ramtek  model  9351  Image  Display  System  selected  for  this 
project  for  the  following  reasons: 

The  modular  architecture  of  the  9351  allowed  the 
selection  of  a  performance  capability  desired  for 
this  application.  Also,  field  expansion  or  modi- 


fication  of  the  system  at  a  later  date  would  be 
easy  to  accomplish. 

-  The  display  refresh  memory  stores  16  bits  per 
picture  element. 

-  The  function  table  provided  for  programmable  defini¬ 
tion  of  output  intensity  and  color. 

The  9351  was  on  GSA  schedule  so  it  could  be  pur¬ 
chased  at  a  substantial  discount. 

Local  service  was  available  from  the  Ramtek  office 
in  Silver  Spring,  Maryland. 

An  existing  software  package  used  for  satellite 
image  processing  was  available  for  the  Ramtek 
which  provided  much  of  the  processing  necessary 
for  display  of  the  microwave  data. 

A  price  quotation  was  obtained  from  Ramtek  to  establish  the 
price  of  the  system.  A  letter  proposal  was  then  submitted  to  the 
sponsor  recommending  Ihc  Ramtek  equipment  and  requesting  addi¬ 
tional  funds  for  its  procurement.  While  our  proposal  was  being 
evaluated,  two  events  occured  which  were  to  impact  this  project. 
First,  Ramtek  instituted  a  10  percent  price  increase  on  all  their 
equipment  and  second,  the  GSA  contract  expired.  Thus,  when  we 
received  authorization  to  proceed  with  acquisition  of  the  equip¬ 
ment,  we  were  facing  a  25  percent  over-run  in  material  cost. 

Over  the  next  few  weeks,  we  proceeded  to  negociate  an  extension 
of  the  original  price  quotation  from  Ramtek  in  exchange  for  an 
extended  delivery  schedule. 

Having  defined  the  image  processing  system,  the  next  task 
was  to  provide  a  hard  copy  capability.  The  most  common  method 
of  making  a  hard  copy  at  this  time  was  to  take  a  photogtaph  of 
the  color  monitor.  However,  this  required  the  services  of  a  photo¬ 
grapher  and  a  darkroom  capable  of  processing  color  prints.  Not 
only  was  this  a  costly  process  but  it  was  not  well  suited  to  a 
research  environment  where  one  is  not  able  to  predict  when  a  hard¬ 
copy  of  the  data  will  be  required.  In  our  search  for  an  alterna¬ 
tive,  we  contacted  the  Polaroid  Corporation  and  were  informed 
that  they  did  not  have  a  camera  which  would  record  the  image 
directly.  However,  we  were  informed  that  two  companies 

Dunn  Instruments,  San  Francisco  CA 

.  Matrix  Instruments,  Milan  MI 

made  comeras  for  use  with  CAT  scanners  and  used  a  newly  developed 
8k  by  11  Polaroid  film.  Subsequently,  we  were  able  to  arrange  for 
a  demo  of  both  camera  systems.  Each  camera  was  fed  an  RGB  video 
signal  from  the  Ramtek  and  a  number  of  photographs  were  taken. 


The  Matrix  camera  seemed  more  reliable  and  produced  a  slightly 
superior  photograph.  The  Matrix  camera  also  had  a  unique  feature 
which  allowed  up  to  16  images  to  be  recorded  on  one  8-2  by  1  i  sheet 
of  film  by  reducing  the  size  of  each  image.  The  capability  of 
producing  multiple  images  on  one  sheet  of  film  was  an  extra  cost 
option.  Each  time  the  number  of  images  was  doubled,  the  price 
increased  by  $2000.00.  In  discussions  with  the  technical  officer, 
it  was  decided  that  two  images  per  sheet  of  film  was  the  optimum 
in  price  versus  capability  consistent  with  their  needs  at  that 
time.  The  camera  could  be  upgraded  later  by  field  modification 
if  the  need  arose. 

On  March  19,  1980,  a  purchase  order  was  issued  to  Matrix 
Corporation  for  a  model  2001  camera  system  for  delivery  in  90 
days.  The  first  delivery  date  was  not  met  and  we  were  informed 
that  production  difficulties  have  caused  the  delivery  schedule 
to  slip.  Many  phone  calls  later,  we  were  informed  that  Matrix 
had  stopped  production  of  the  2001  camera  and  was  in  the  process 
of  redesigning  the  entire  system.  The  new  camera,  called  a  model 
4007,  was  to  incorporate  all  the  optional  features  of  the  2000 
series  camera  plus  have  increased  reliability  and  ease  of  opera¬ 
tion.  If  we  would  agree  to  a  six  month  delivery,  Matrix  proposed 
to  deliver  the  new  model  for  the  same  price  as  the  model  2001. 

In  discussions  with  the  technical  officer,  it  was  decided  that 
the  six  month  delay  was  accptable  in  return  for  a  camera  system 
with  much  improved  performance.  As  a  result  of  this  change,  the 
government  received  a  camera  system  worth  $4000  mote  that  the 
price  they  were  charged.  Unfortunately,  the  schedule  on  the  new 
camera  also  slipped  and  it  was  not  delivered  until  April  3,  1981. 
It  should  be  noted  that  the  Matrix  camera  has  lived  up  to  all 
the  expectations  and  is  considered  well  worth  the  extended  time 
necessary  to  obtain  it. 

B .  System  Description 

Figure  1  is  a  detailed  block  diagram  of  the  RM-9351  System. 
The  various  subsystems  are  briefly  described  in  the  following 
paragraphs . 

Host  Parallel  Link  -  The  host  parallel  link  provides  the  high 
speed  (up  to  600  RHz)  bidirectional,  16-bit  parallel  communica¬ 
tions  between  the  host  computer  and  the  display  system.  Four 
external  interrupts  are  implemented.  The  computer  uses  direct 
memory  access  (DMA)  to  communicate  with  the  Ramtek  GPIF,  which 
is  typically  contained  within  the  control /video  board. 

Internal  Processor  Bus  -  The  internal  processor  bus  connects 
the  host  parallel  link,  Z80  microprocessor,  registers,  memories, 
video  generators,  and  interfaces  for  optional  peripheral  equip¬ 
ment.  This  bus  provides  high-speed  device-to-device  communica¬ 
tion  and  implements  nonprocessor  DMA. 

Display  Address  and  Data  Registers  -  Display  address  and  data 
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registers  connect  the  bus  and  refresh  memory  for  image  generation 
and  retrieval.  This  communication  is  in  the  form  of  a  16-bit 
word  per  pixel,  with  up  to  16  bits  being  written  in  the  z-axis. 

The  X-  and  Y-  addresses  of  the  pixels  are  incremented  as  prescribed 
by  the  control  registers  and  logic  in  the  display  generator.  The 
display  processor  interface  is  optimized  with  respect  to  its  in¬ 
ternal  algorithms  for  generating  character  and  graphics  data. 

Graphic  Display  System  -  Graphics  data  is  written  under  the 
direction  or  the  current  operating  point  (X-Y  address),  that  is, 
up  or  down  and/or  left  ot  right.  The  color  or  intensity  of  font, 
raster,  and  graphics  data  (Z-axis)  is  assigned  by  the  Z80  micro¬ 
processor.  Foreground,  background,  reversal,  and  writing  mode 
(replacement  versus  additive)  may  be  specified.  The  Z80  micro¬ 
processor  is  a  powerful  tool  for  implementing  not  only  these  modes 
of  writing  but  any  other  application-dependent  mode. 

Video  Monitor  -  The  video  monitor  decodes  the  generated  video 
s igna 1  and  displays  Lho  image  by  driving  one  or  more  cathode-ray 
beams  in  raster  fashion.  Tube  refresh  time  is  at  30~Hz  for  momes- 
tic  systems. 

Display  Processor  -  This  display  processor  interprets  display- 
instruction  information  and  presides  over  the  bus.  Secondary 
functions  include  character  generation,  vector  generation,  plot 
generation,  raster  mode,  and  raster  margining.  A  Z80  micropro¬ 
cessor  with  1,024  bytes  random-acccss-mcmory  (RAM)  and  5,120 
bytes  of  PROM  is  included.  A  basic  instruction  set  is  implemented, 
which  provides  imaging,  graphics,  and  text-generating  functions. 

Refresh  Memory  -  The  refresh  memory  provides  sixteen  bits  of 
storage  for  refreshing  each  pixel  on  the  CRT.  The  refresh  memory 
contains  two  memory  boards,  each  board  containing  eight  sections 
of  memory.  Each  section  contains  sixteen  16K  MOS  RAMs  that  offer 
512-line  by  512-element  refresh  capability  in  the  RM-9351  model. 

One  plane  of  memory  stores  one  bit  of  storage  per  pixel.  For  ex¬ 
ample,  a  refresh  memory  board  in  the  RM-9351  contains  sixteen 
sections  of  memory,  i.e.,  one  plane  per  section.  Thus,  the  RM- 
9351  offers  16  planes  or  16  bits  of  storage  per  pixel. 

Video  Generation  -  The  video  generation  section  of  the  control 
board  can  be  modified  through  PROM  coding  to  yield  a  multitude 
of  possible  configurations.  Figure  2  shows  the  video  generation 
network  contained  on  the  control /video  board.  By  utilizing  only 
the  basic  control/video  board  and  one  eight-plane  memory  board, 
the  user  can  specify  the  PROM  code  so  that  the  basic  four  video 
amps  provide  him  with  any  of  the  following  systems: 

a.  Four  planes  (plane  0  through  plane  3)  patched  to 
four  video  amps  to  yield  four  black-and-white  dis¬ 
plays.  This  is  the  standard  configuration. 

b.  Six  planes  patched  to  three  video  amps  (plane  0 
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through  plane  5)  to  yield  64  colors  to  drive  an 
RGB  monitor.  The  remaining  video  can  be  used  to 
drive  a  composite  signal  of  planes  0  through  5  to 
a  B-W  monitor  which  can  provide  four  grayscale 
levels . 

C.  Eight  planes  patched  to  provide  64  colors  on  an 
RGB  monitor  (planes  0  through  5)  with  two  planes 
provided  to  yield  two  overlays  (planes  6  and  7). 

The  remaining  video  amp  may  drive,  using  planes 
0  through  7,  a  composite  black-and-white  picture. 

Enhancement  tables  are  video  lookup  tables  that  allow  inter¬ 
active  pseudocolor  or  gray-scale  translation  of  images  stored  in 
refresh  memory.  True  color  results  when  a  color  monitor  is  driven 
directly  by  memory  data;  pseudocolor  results  when  memory  data 
is  manipulated  with  a  RAM  before  being  sent  to  a  color  monitor. 

When  equiped  with  a  lookup  table,  the  stored  refresh  data  is 
treated  as  an  address  to  the  lookup  table,  whiph  is  host  program¬ 
mable;  that  is,  as  each  pixel  is  scanned  from  the  refresh  memory 
for  video  presentation,  the  contents  of  the  corresponding  cell  is 
retrieved  from  the  lookup  table  and  this  data  is  passed  to  the 
digital-to-analog  converters  and  video  amplifiers,  instead  of  the 
refresh  data  itself.  Thus,  the  refresh  data  addresses  a  host  pro¬ 
grammable  lookup  table  that  assigns  the  output  intensity  or  color. 

Lookup  tables  are  most  often  used  in  imaging  applications. 

Figure  3  illustrates  a  commonly  used  enhancement  procedure  termed 
windowing,  or  density  (level)  slicing.  Here,  a  specified  range 
or  group  of  contiguous  image  intensity  values  are  fitted  to  the 
available  spectrum  of  output  intensity  levels.  Pixels  having 
values  beneath  or  above  the  specified  range  (or  window)  are  trans¬ 
lated  to  black  or  white,  respectively;  while  pixels  within  the 
window  are  translated  to  an  appropriate  gray  level.  Thus,  the 
observer's  attention  is  focused  upon  the  window,  and  he  can  more 
easily  distinguish  between  what  were  minor  or  negligible  inten¬ 
sity  differences  in  the  original  image.  More  important,  the 
user  is  able  to  manipulate  this  window  interactively  by  reloading 
the  lookup  table,  and  without  affecting  the  image  in  refresh 
memory.  The  integrity  of  the  original  image  is  thus  maintained 
while  the  visual  presentation  is  varied  to  suite  the  needs  of  the 
observer. 

Figure  3  illustrates  the  effect  of  a  relatively  simple  con¬ 
trast  enhancement  algorithm.  It  is  important  to  note  that  the 
equipment  supports  sophisticated  algorithms.  For  example,  a 
gamma  corrected  output  can  be  achieved  by  loading  a  non-linear 
function  into  the  lookup  table,  or  a  pre-defined  or  computed  set 
of  pseudocolors  might  be  assigned  to  a  gray  scale  image,  etc. 

The  enhancement  table  consists  of  2,048  cells  of  nine  or 
thirteen  bits.  Thus,  up  to  11  bits  of  image  data  (in  refresh  memory) 
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can  address  the  table.  For  gray  scale  applications,  bit  2q 
through  2y  of  the  cells  are  output  to  a  single  8-d igi tal-to-ana lo 
converter  (DAC) .  Thus  any  of  256  gray  level  intensities  (i)  may 
be  produced,  i.e.,  1  =  N/255  where  0  <  N  ^  255.  For  color  appli 

cations ,  bits  20  through  23,  24,  through  2  7>  and  2g  through  2n 
are  output  to  three  sepcrate  4-bit  DAC ' s  corresponding  to  the 
red  (R),  green  (G) ,  and  blue  (B)  primary  inputs  to  an  RGB  color 
monitor.  Thus,  any  of  4,096  colors  (C)  can  be  produced,  i.e., 

C  =  2g  x  R  +  24  x  G  +  2g  x  b  where  0  <  R,G,B  15.  Whether  gray 
scale  color,  bit  2^,  when  non-zero,  causes  the  corresponding 
color  or  intensity  to  blink  at  a  2  Hz  frequency. 

C .  Image  System  Operation 

The  operation  of  the  image  system  was  designed  to  provide 
the  user  with  a  very  capable  system  that  is  simple  to  operate. 

All  communication  between  the  system  and  user  takes  place  through 
the  computer  CRT  console.  A  menu  of  the  system  commands  is  dis¬ 
played  on  the  CRT  as  an  aid  to  the  operator.  To  select  a  parti¬ 
cular  function,  the  operator  need  only  type  the  first  two  charac¬ 
ters  of  the  particular  command.  The  image  system  commands  are 
listed  below  along  with  a  brief  description  of  their  function. 

STOP,  END,  OR  EXIT 

terminates  program  execution. 

INITIALIZE 


Program  is  initialized.  All  paramters  are  set  to  their  default 
values.  Next  input  command  may  now  be  chosen  from  list  on  CRT. 

GRID 


This  gives  the  user  the  capability  of  varying  either  the  picture 
or  the  grid  intensity.  The  program  will  prompt  the  user  with 
Input  grid  in  (0  to  1).  The  user  must  input  the  desired  grid 
brightness  by  punching  in  a  number  between  0.0  (black)  and  1.0 
(white).  A  number  greater  than  1.0  will  cause  the  grid  to  become 
white  and  the  brightness  of  the  picture  will  decrease  by  an  in¬ 
versely  proportionate  amount.  I.e.,  if  a  value  of  two  is  input 
the  picture  will  become  half  as  bright  while  the  grid  will  become 
fully  white.  If  a  value  of  .5  is  input  the  picture  will  be  at 
maximum  brightness  while  the  grid  will  become  half  as  bright. 

To  return  to  R0L0R  any  of  the  input  commands  on  the  CRT  may  be 
selected . 

TRANSFER 


This  command  transfers  a  new  picture  from  the  disc  to  the  RAMTEK 
display.  The  program  will  prompt  the  user  with  Picture  disc 
Logical  Unit.  Punch  in  the  LU  number  and  press  RETURN  key.  The 
program  will  then  prompt  "Picture  number".  Input  the  number  of 


the  picture  you  want  transferred  and  press  RETURN  key.  The  appro- 
riate  picture  will  be  sent  to  display  system.  On  the  CRT  the 
prompt  picture  nummber  is  returned  in  case  another  picture  trans¬ 
fer  is  required.  Any  ROLOR's  commands  may  also  be  input  to  return 
control  to  ROLOR  and  then  to  go  to  the  specified  command. 

COLOR 


This  allows  the  user  to  select  different  color  schemes  for  the 
pictures.  The  choice  is  varied  and  listed  below.  The  program  will 
prompt  the  user  with  choice  (A  through  K,  BW,  A1  A9  or  B1  B9) . 
Choose  a  color  scale  (from  the  list  below)  and  input  your  choice. 
The  color  schme  of  the  picture  will  change  at  once.  The  input 
prompt  will  remain  in  case  another  color  change  is  required.  Any 
of  ROLOR's  commands  may  also  be  input  to  return  control  to  ROLOR 
and  then  go  to  the  specified  command. 

Key  to  color  codes: 

A-  Spectral  (ROYGBIV) 

B-  Half  contouring  of  A 
C-  Full  contouring  of  A 

D-  Green,  Blue,  Black,  Green,  Yellow,  Red,  Pink 

E-  From  Bull  of  Amer  Meteorlogical  Soc .  Vol.  52,  #9  (Sept  71) 

F-  For  NIMBUS  5  -  ESMR  pictures 
G-  Yellow,  Dark  Orange,  Dark  Blue,  Light  Blue 
H-  Blue,  Green,  Violet,  Dark  Orange  White 
I-  Blue,  Green,  Yellow,  Red,  Pink,  White 
J-  Variation  of  I 

K-  Pink,  Orange,  Yellow,  Green,  Blue 
L-  Microwave  pseudocolor  scale 
BW-  Black  and  white 

A1  through  A9-  Contour  intervals  in  green 

A1  contours  one  interval  (center  of  scale) 

A2  contours  two  equally  spaced  intervals,  etc. 

Position  of  intervals  contoured  can  be  shifted  by 
changing  the  color  scale  limits  using  the  SC  command 

B1  through  B9-  Same  as  A1  through  A9 ,  except  black  and  white 

BACKGROUND 


This  changes  the  background  color  for  the  title  and  color  bar. 
The  background  color  is  given  by  three  values,  each  ranging  from 
0  to  1,  which  specify  the  intensities  of  red,  green  and  blue 
respectively  in  the  background. 

e.g.  0,0,0  for  black  background 

0,1,0  for  green  background 

.4,0,. 6  for  purple  background 

when  the  program  prompts  Input  red,  green,  blue  (0  tol)  values 
then  input  your  selection  and  press  the  RETURN  key.  The  color 
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change  is  immediate  on  the  RAMTEK  display  system.  The  prompt 
will  reappear  in  case  another  color  change  is  required.  If  not, 
any  of  ROLOR's  input  commands  (listed  on  CUT)  may  be  used. 
Control  is  returned  to  ROLOR  which  then  executes  the  specified 
command . 


BLOW-UP 

This  allows  the  user  to  enlarge  a  portion  of  the  picture.  The 
program  will  prompt  the  operator  with  "Picture  disc  Logical  Unit." 
Input  LU  number  of  the  picture  disc.  Press  the  RETURN  key.  The 
program  will  then  prompt  "Input  line  #,  column  #,  factor,  (Piet  #).' 
The  line  #  and  the  column  #  specify  the  center  of  the  area  that 
is  to  be  enlarged.  Values  from  1  to  511  may  be  input  for  the  two 
parameters.  "Factor"  is  the  amount  by  which  the  picture  is  to  be 
enlarged.  I.e.,  a  factor  of  two  will  make  the  picture  twice  as 
large.  The  picture  #  is  an  optional  parameter  for  the  number  of 
picture  that  is  to  be  enlarged.  If  not  specified  it  defaults  to 
the  value  of  the  current  picture.  Input  your  selection  and  press 
RETURN  key.  The  prompt  specified  above  will  reappear  in  case  an¬ 
other  change  is  needed.  If  not,  any  of  ROLOR's  commands  may  also 
be  used.  That  command  will  then  be  executed. 


DATA  or  SCALE 

This  rpvides  the  user  with  the  ability  to  change  the  limits  of 
the  data  scale  and  has  a  kaleidoscope  and  upkill  facility.  The 
user  must  supply  four  parameters,  spearated  by  commas,  when  the 
program  gives  the  prompt  "Input  scaling  for  white  1/black  =  0, 
upkill,  klido" 


1st  para  -  Upper  limit  of  data  scale.  If  less  than  1  the  data 

with  values  beyond  the  limit  will  be  displayed  as  white 
(or  the  current  background  color). 

2nd  para  -  Lower  limit  of  data  scale.  If  greater  than  0  the  data 

with  values  below  the  limit  will  be  displayed  as  black. 

3rd  para  -  Upkill  (0  or  1).  If  the  value  is  1  data  with  values 

above  upper  limit  is  displayed  as  black.  If  the  value 

is  0  the  data  is  white.  It  is  an  optional  parameter 
and  defaults  to  0. 


4th  para  -  Kaleidoscope  facility.  This  causes  the  picture  to 
cycle  through  brightness  scales  automatically.  The 
larger  the  number  the  slower  the  cycling.  The  length 
of  time  for  which  it  runs  is  fixed  and  operator  must 
wait  for  the  whole  sequence  to  end  before  issuing  an¬ 
other  command. 


Default  values  are  1,0,0,  no  kaleidoscope.  The  third  para¬ 
meter  is  optional.  If  not  specified  it  equals  0.  The  fourth 
parameter  is  also  optional.  If  not  specified  there  is  no 
kaleidoscope . 
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If  the  lower  limit  is  greater  than  the  upper  limit,  the  color 
scales  are  complimented. 

If  an  error  occurs  control  is  returned  to  ROLOR  and  you  must 
specify  the  next  input  command  from  those  listed  on  CRT. 

ARCHIVE 

This  command  references  another  program,  program  RATPT.  It  pro¬ 
vides  the  user  with  the  ability  of  transferring  a  complete  512  x 
512  picture  between  the  disc,  the  mag  tape  or  the  RAMTEK.  To 
initiate  a  transfer,  the  operator  must  select  one  of  the  command 
options  displayed  on  the  CRT,  type  in  the  two  letter  mnemonic 
for  the  command,  and  press  the  RETURN  key. 

OVERLAY 

This  command  references  another  program,  program  OVTPR.  It  pro¬ 
vides  the  user  with  the  capability  of  changing  the  grid  and  back¬ 
ground  colors  of  a  picture  to  suit  the  picture  requirements. 


SECTION  III 


IMAGE  PROCESSING  PROGRAMS 
A.  CUBIC  SPLINE  INTERPOLATION 


The  microwave  scanning  system  produces  a  64  x  64  pixel  image. 
The  RAMTEK  image  system  supports  a  512  x  512  pixel  image,  so  the 
raw  image  would  occupy  only  one  eighth  of  the  screen,  and  little 
detail  would  be  visible.  There  are  numerous  ways  to  enlarge  an 
image,  such  as  simple  pixel  replication  or  bilinear  interpola- 
tion.  While  these  methods  are  easier  to  implement  than  cubic 
spline  interpolation,  they  have  the  disadvantage  of  degrading 
fine  detail  in  the  image.  Thus,  it  is  worth  the  effort  to  inter¬ 
polate  with  cubic  spline  functions. 

The  particular  method  of  cubic  spline  interpolation  which 
was  finally  selected  for  implementation  is  outlined  in  Numerical 
Methods  by  Robert  Hornback.  This  source  describes  how  to  fit  a 
cubic  spline  to  a  set  of  points  (xi,  f(xi>).  This  is  done  as 
follows.  Between  each  pair  of  adjacent  points  xi  and  x£+i  it  is 
necessary  to  find  a  cubic  polynomial  which  passes  through  (xi, 
f(xi))  and  (xi  +  i , f (x£  +  i ) )  .  This  polynomial  is  denoted  by  Fi(x)  = 
ao+a].x+a2x2+a3x3  for  xi  <  x  <  xi+i.  There  are  4  unknown  constants 
^i(xi+i)  =  f  (x£+i).  The  remaining  constraints  are  imposed  by 
requiring  that  the  first  and  second  derivatives  of  F£  match  those 
of  the  polynomial  Fi_i  used  on  the  previous  interval.  This  gives 
the  cubic  spline  interpolator  its  characteristic  smoothness.  To 
actually  calculate  the  cubics  Fi,  it  is  first  necessary  to  calcu¬ 
late  the  second  derivatives  g”^)  which  can  be  found  by  solving 
the  following  set  of  simultaneous  equations: 


A  Xi-x 

g"(xi-i)  + 

2(xi+1  _  Xi_! ) 

_  A  xi 

A  xi 

g"(xi)  +  g'^Xi-i) 


+  6 


f(xi+l)  -  f(xi)  f(Xi>  _  fCx^) 


(  &  xi)2 

where  x£  =  xi+1  _  x£. 


(  Axi)  (  Axi-i) 


(i  =  1,2,.. ,n-l) 


Note  that  this  procedure  only  n-1  equations  in  n  +  1  unknowns 
g " ( xq ) . .  g"(xn).  The  2  additional  equations  are  obtained  by 
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specifying  conditions  on  g"(x0)  and  g"xn) .  For  this  purpose, 
take  g"(x0)  =  0  and  g"(xn)  =0. 

For  image  processing,  consider  the  x£  to  be  at  integer  values 
and  the  functional  values  to  be  the  pixel  intensities.  Thus  x^  =1, 
and  the  equations  are  simplified.  The  functions  Fj  can  now  be 
written  in  terms  of  g"  and  f: 


Fi(*>  » 

6 


(xi  +  l  x)3  _  (xi  +  i  x}]  +  g"(x£+j)  x  _  xj)3  —  (x~x^ ] 


+  fixi)  (xi+1-  x)  +  f(xi+l)  (x  _  x.; 


As  an  example 
sional  edge  which 
512,512.  In  this 

512. 


,  consider  the  effect  of  expanding  a  one  dimen- 
consists  of  five  pixels  with  intensities  000 
case,  x0  =  0,...,x4  .  4  and  f(x0)  =0,...,  f^)  m 


We  get  the  following  conditions  on  t>"  • 

1  II  I  .  °  * 


hence  g"(l) 

g"(2)  =i 

g"(3)  56 


g"(l) 

g"(2) 

0 

=6  512 

g  (3) 

-512 

-5k 

12k 

-19k 


512  I  ,  g"(0),  g" (4)  =  0 

5k 

2k 

ak  where  k  =  512 


1 

and  the  equations  are: 

F0(x)  =  — ^jr  (x^  -  x) 

F1(X)  =  "H  f2  -  x)3  -  (2  -  x)]  = 


H 


.[(X  "  l)3  -  (x  -  1)  ] 


?2(X)  =  ^  [(3_X)3-  (3-X3  -^[x-2)3_(x.2)]  + 


k(x  -  2) 


F3(x) - 19k  [(4  -  x)3  _  (4  -  X)  ^  +  k 

To  interpolate  one  pixel  in  between  each  original  pixel  we 
must  compute  F0(.5J,  Fl(l.5>,  ^(2.5),  and  F3(3 Jj. 


.5)  =  .033k  =  17 

f  1  ( 1 . 5 )  =  -.047k  =  -24  =  0  (lor  image  processing) 

F2(2 . 5)  =  .547  =  280 
F3(3.5)  =  1.127k  =  577 

The  expanded  edge  is  slightly  enhanced  and  smoothed,  as 
shown  in  Figure  4. 

Since  the  cubic  spline  interpolator  applys  to  one  dimensional 
arrays  of  data  and  an  image  is  two-dimensional,  it  is  essentially 
necessary  to  make  two  passes  with  the  interpolator.  First,  treat 
each  column  of  the  image  as  a  set  of  data  and  interpolate  the  re¬ 
quired  number  of  pixels  in  each  column.  Then,  treat  each  original 
plus  interpolated  row  of  the  image  as  a  one  dimensional  data  set 
and  interpolate  horizontally. 

Conceptually,  the  mthod  is  shown  in  Figure  5. 

In  practice,  to  conserve  memory  requirements  not  all  the 
vertical  interpolation  is  done  at  once.  Instead,  all  the  infor¬ 
mation  needed  to  construct  the  interpolating  polynomials  for  each 
column  is  computed  and  stored.  Then  vertical  interpolation  is 
performed  on  two  adjacent  column  pixels  for  every  column,  followed 
by  horizontal  interpolation  across  the  original  and  interpolated 
rows.  Then  the  data  is  flushed  to  the  Ramtek  and  work  begins 
on  the  next  adjacent  column  pixels. 

To  allow  the  user  maximum  flexibility,  this  implementation 
of  the  cubic  spline  interpolation  is  not  hard  coded  to  any  parti¬ 
cular  expansion  factor.  Due  to  screen  size  limitations  it  is  not 
proactical  to  insert  more  than  7  pixels  between  the  originals, 
but  the  program  can  handle  any  expansion  factor  between  1  and  8, 
inclusive . 


B.  TWO  SEGMENT  ENHANCEMENT  CURVE 

The  addition  of  a  two  segment  enhancement  curve  to  the  RAMTEK 
9351  system  was  accomplished  by  implementing  a  new  subroutine  and 
operator  calls  to  the  previously  existing  RAMTEK  software  package. 
The  implementation  allows  a  choice  of  single  line  or  dual  line 
enhancement  curves  by  appending  additional  information  to  the 
enhancement  operator  query.  Because  the  enhancement  curve  must 
be  single  valued,  various  tests  are  included  to  prevent  the  gen¬ 
eration  of  meaningless  curves.  Internally,  the  software  multi¬ 
plexes  the  requested  color  scale  and  the  enhancement  curve  to 
produce  the  video  lookup  table  bit  format,  which  is  then  trans¬ 
ferred  to  the  Image  System.  The  enhancement  curve  with  descrip¬ 
tions  is  shown  in  Figure  6. 
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0001  FTN4  ,L  APPENDIX 

00  02  c*** **************************************** ********** 


0003 

C* 

Uritten  by  Ginny  Kalb 

* 

0004 

c* 

Technology  U.  S.  A. 

* 

0005 

c* 

P.  0.  Box  55333  Ph  <301)  292-2592 

* 

00  06 

c* 

Oxon  Hill  Station 

* 

0007 

c* 

Uashington  D.  C.  20022 

* 

0008 

c* 

Rev .  29  Jan  1981 

* 

0009  C*** **  ************************************************************  ** 

0010  C  This  program  uniformly  expands  a  square  picture  in  both  dim- 
0011  C  ensions  by  using  a  cubic  spline  i nterpo 1  at i on .  This  is  done 
0012  C  by  first  interpolating  along  the  columns  of  the  input  picture, 

0013  C  then  across  the  original  plus  interpolated  rows.  The  resultant 
0014  C  image  is  displayed  on  the  RAMTEK  or  COMTAL  screen  and  optionally 
0015  C  stored  on  disc.  Padding  with'zeros  is  used  to  obtain  a  full 
0016  C  screen  output  picture. 

0017  C 

0018  C  NOTE:  Several  array  sizes  will  need  to  be  altered  for  differ- 
0019  C  ent  input  picture  dimensions,  as  well  as  parameter  IDIM  = 

0020  C  pixels/side  of  input  picture.  The  array  sizes  are  a  function 
0021  C  of  IDIM  and  are  located  in  COMMON  blocks  INPUT,  UNIQ,  and  RROW. 
0022  C  Subroutine  TREAD  is  of  course  input-unique  and  must  be  altered 
0023  C  for  a  new  input  picture  size. 

0024  C  Array  sizes  in  subroutine  INVRS,  used  by  SOLVE,  must  be  changed 
0025  C  to  A< IDIM-2, IDIM-2)  and  J<IDIM+19>. 

0026  C 

0027  C  COMMON  /'INPUT/' IBUF<  IDIM,  IDIM  ) 

0028  C  COMMON  /UNIQ/NPTS ,  IDIM,  CCO£F<  IDIM,  4  >,  DELTA,  ISHF,  XSTEP , 

0029  C  DERIV< IDIM, IDIM-2>,T< IDIM-2, IDIM-2) 

0030  C  COMMON  /RR0U/'A<  4  >,  S<  IDIM-2  ),  Y<  IDIM-2  > 

0031  C 

0032  C  The  number  of  interpolated  points  is  a  user  option  at  run  time. 
0033  C  The  default  is  3  which  expands  a  64x64  picture  to  253x253.  The 
0034  C  maximum  value  for  a  64x64  input  picture  is  7  because  after  that 
0035  C  the  output  picture  would  exceed  512x512. 

0  036  C  '•¥  +  *  ’***  *  *  *1**4*  *  *  *  *•*  *4*  %  *  ♦  sfnfc  *  *  *  *  4c*  ***%%*  *  %  :*  *  *  *  *  %  *  *  sfr  %  *  ■*  *  *  *-  *  *  *  *  *  ♦  * 

0037  C^Y 
0038  C 

0039  PROGRAM  EXPND 

0040  COMMON  / INPUT/ IBUF< 64 , 64  ) 

0041  COMMON  /0UTPT/I0PT,LUD, IPIC,L, JBUF< 5 1 4  ) , NPRT, LU3, I TK , ISECT 

0042  COMMON  /UN  I Q/NPTS , I D IM , CC0EF< 64 , 4 > , DELTA , I  SHF, XSTEP , 

0043  *  DERIV<64,62),T<62,62) 

0044  COMMON  /RR0U/A< 4 >, S< 62 >, Y< 62 > 

0045  DIMENSION  KBUF<512> 

0046  EQUIVALENCE  < KBUF< 1  ), J8UF< 3 > > 

0047  DIMENSION  IP<5> 

0048  C 

0049  CALL  RMPAR< IP  ) 

0050  URITE< IP, 1  ) 

0051  1  F0RMAT<  "tm^htJ",  "£S,a  1  7r  OC^&dOEXPANS  I  ON  VIA  CUBIC  SPLINE") 

0052  UR  I TE< IP , 2  ) 

0053  2  F0RMAT<  "  fcta  1  9r  OC^KfcU.dCNOTE^S.dl?  program  defaults  to  RAMTEK", 

0054  *  "  screen  on  LU  1  2  ,  " ,  "^Sca2  Or  OC^KDo  you  want  to  change  either", 

0055  *  "  default?  Y/N  ) 

0056  READ< IP, 30  >  ITEMP 

0057  IF< ITEMP. EQ. 1HN>  GO  TO  6 

0058  UR  I  TE<  I P  ,  3  ) 


0059 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076 
0077 
0078 
0079 
0080 
0081 
0082 
0083 
0034 
0085 
0086 
0087 
0088 
0089 
0090 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
0098 
0099 
•  0100 
0101 
01  02 
01  03 
01  04 
01  05 
01  06 
0107 
0108 
01  09 
01  10 
01  1  1 
01  12 
0113 
0114 
01  15 
01  16 
01  17 
01  18 


3  FORMAK  "5;8<a20r0C$:K"  ,  “t*«a1  9rOC^fcdC  .  t**d8input  LidCCO", 

*  "tide  for  COMTAL  or  t*«dCRAt&de  for  RAHTEK  > 

READ<IP,4>  NPRT 

4  FORMAK  A2) 

URI TE< IP ,  5  > 

5  FORMAK  utia20r0CtK“,  "tia19rOCtKt*.dC.tid0inpufc  LU  #  _'*  > 

READ< IP,*)  LU3 

6. CONTINUE 
UIRITE<  IP, 7) 

IJRITE<  IP,  8  > 

7  FORMAK  Mtia19r0CtKtidC.tidespecify  expansion  factoi - ", 

*  “  e.g.  2  will  double  input  picture  which  is  64x64  ") 

8  FORMAK  "t««a20r0CtK<or  just  RETURN  and  use  the  default  of  4  >"  > 
I TEMP  ■  -1 

REA0< IP, * )  ITEMP 

IF< ITEMP.GE. 0)  NPTS  -  ITEMP-1 

IF<<< IOIM-1  )*NPTS+ IDIM ) . GT . 5 1 2 ) . GO  TO  995 

URITE< IP, 1 0 ) 

10  FORMAK  "II^OrOCtK",  "fci,a19rOCtK£6.dC.fc8,d0specify  input  ", 

*  "picture  tape  unit  _M ) 

REAO< IP, *  >  MT 

C 

CALL  INVRS<  T , 62, OTNRM, DETM ) 

C 

15  URITE< IP, 20) 

20  FORMATC “tia20r0CtK", utia19r0CtKtidC.tid8do  you  want  ", 

*  “to  save  output  picture?  Y/N  > 

READ< IP, 30  >  IOPT 

30  FORMAT<  A 1  ) 

IF< IOPT.EQ. 1HN>  GO  TO  50 
URITE< IP, 40) 

40  FORMAK "tia19rOCtKtidC.tid8specify  output  disc  LU  and", 

*  "  picture  number  “) 

READ< IP, *  )  LUO, IPIC 

C 

C  Verify  that  the  .specif ied  .output  disc  is  legitimate 

C 

CALL  .CHECK< IP,LU D) 

C 

50  CONTINUE 
C 

C  Initialization  of  parameters 
C 

L  -  -1 

ISHF  -  <  512-<< IDIM-1 >*NPTS+IDIM>>^2 
XSTEP  -  <  NPTS+1 .0  >*DELTA 
ICNUD  -  MT+100B 
C 

C  If  using  RAMTEK,  issue.a  reset  command 
C 

IF<  NPRT . EQ , 2HCO  >  GO  TO  65 
CALL  EXEC< 2, LU3, 2400B,  1  ) 

C 

C  Output  leading  rows  of  zero 
C 

65  LIM*<512-<<  IOIM-1  >*HPTS+  IDIM  >  V2 
LOSK  »  <480-<<  IDIM-1  )*NPTS  +  IDIM)V2 
IF< NPRT .  EQ . 2HC0 )  LIM  -  LDSK 
IF<LIM . LE . 0 >  GO  TO  72 
DO  70  K-1,LIM 
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Of  19 

IF<k\LE.LDSK)  L  »  K-1 

0120 

CALL  DSPLAY 

0121 

IF< I0PT.EQ, 1HN)  GO  TO  70 

0122 

IF< K . GT . LDSK >  GO  TO  70 

0123 

CALL  EXEC<2,LUD,KBUF,512, ITK,ISECT> 

0 1  24 

70 

CONTINUE 

0125 

C 

0126 

C  Read  in  entire  picture 

0127 

C 

0128 

72 

CONTINUE 

0129 

00  75  1*1 , I DIM 

0130 

CALL  TREA0< ICNWD, I, IP > 

0131 

75 

CONTINUE 

0132 

C 

0133 

C. Solve  for  column-determined  second  derivatives 

0134 

C 

0135 

DO  100  J-1 , IDIM 

0136 

CALL  S0LV£<  0, J  > 

0137 

100 

CONTINUE 

0138 

C 

0139 

C  Interpolate  across  first  row 

0140 

C 

0141 

DO  21 0  J*1 , IDIM 

0142 

INDEX  «  ISHF+<NPTS+1 >+J-NPTS 

0143 

KBUF< INDEX  >  -  IBUF<f,J> 

0144 

210 

CONTINUE 

0145 

CALL  ROUk  1  > 

0146 

C 

0147 

C  Loop 

on  remaining  rows,  doing  column  followed  by 

0148 

C  i  nte 

rpolation 

0149 

C 

0150 

DO  350  1*2, IDIM 

0151 

N  =  1-1 

0152 

DO  31 0  J-1 , IDIM 

0153 

CALL  GETCO<  0,  J,N> 

0154 

31  0 

CONTINUE 

0155 

C 

0156 

IF< NPTS . EQ . 0 )  GO  TO  331 

0157 

DO  330  K* 1 , NPTS 

0158 

X  «  <N-1 )*XSTEP+K*DELTA 

0159 

DO  320  J-1, IDIM 

0160 

INDEX  -  ISHF+<NPTS+1 >*J-NPTS 

0161 

A< 1  )  »  CC0EF< J, 1 > 

0162 

A<2>  »  CC0EF<  J,  2  > 

0163 

A<3)  *  CC0EF<  J,  3  ) 

0164 

A<4)  -  CCOEF<  J,  4  ) 

0165 

K8UF< INDEX)  ■  F<A,X,N> 

0166 

IBUF< 1-1 , J  >  »  KBUF< INDEX) 

0167 

320 

CONTINUE 

Ot  68 

CALL  R0U<  N  > 

0169 

330 

CONTINUE 

0170 

331 

DO  340  J-1 , IDIM 

0171 

INDEX  -  ISHF+<NPTS+1 >*J-NPTS 

0172 

KBUF< INDEX)  -  IBUF<I,J) 

0173 

34  0 

CONTINUE 

0174 

CALL  R0U< I > 

0175 

350 

CONTINUE 

0176 

C 

0177 

C  Write  trailing  rows  of  zero 

0178 

C 
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0179 
0180 
0181 
0182 
0183 
0184 
0185 
0186 
0187 
0188 
0189 
0190 
0191 
0192 
0193 
0194 
0195 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
02  04 
0205 
0206 
0207 


0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 


DO  400  1*1,512 
K8UF< I >  -  0 
400  CONTINUE 
LSAVE  «  L+l 
DO  410  L»LSAVE,51 1 
CALL  DSPLAY 

IF< IOPT.EQ. INN)  GO  TO  410 
IFCL.GT,479)  GO  TO  410 
CALL  EX£C<  2, LUD, KBUF, 512, I TK, I  SECT) 

410  CONTINUE 

URITE< IP, 1 ) 

WRITE< IP, 450  > 

450  FORMAT< "fcta 1 9r OCtKi t  is  finished") 

WRITE< IP, 455) 

455  FORMAT< M^ia20r OC^K^&dC . ^&d8do  you  want  to  expand  the 

*  "next  picture?  Y/'N  > 

READ< IP , 30  >  ITEMP 
IF< ITEMP. EQ , 1 HN  >  STOP 
ICNUD  -  MT+1300B 
CALL  EXEC<3, ICNWD) 

CO  TO  15 

995  URITE< IP, 996) 

996  FORMATC "(tal 9r OC^Kinval id  value  for  expansion  factor" 

*  "(tia20r OC^&KTry  a  value  between  1  and  8,  inclusive"  > 

STOP 
END 

SUBROUTINE  R0W< I  ) 

This  subroutine  interpolates  NPTS  points  between  the  gray  levels 
in  row. I  of  the  input. buffer  2BUF.  This  is  either  the  original 
row  or  an  interpolated-between-columns  row. 

The  resultant  row  is  then  displayed  and  optionally  stored  on  disc 

COMMON  /'OUTPT/IOPT,  LUD,  IPIC,L,  JBUF<  514),  NPRT ,  LU3,  ITK,  ISECT 
COMMON  7UNIQ7NPTS, IDIM, CCOEFC  64, 4  >, DELTA, ISHF, XSTEP, 

*  DERIV<64,62),T<62,62) 

COMMON  /RR0U/'A<  4>,S<62),Y<62> 

DIMENSION  KBUF<  512) 

EQUIVALENCE  < KBUF< 1 ), JBUF<  3 ) > 

IF< NPTS . EG . 0  )  GO  TO  .201 
CALL  SOLVE< 1,0) 

LIM  »  IDIM-1 
DO  200  J»1  AIM- 
CALL  GETCO< I, 0, J> 

DO  1 00  K*  1  , NPTS 

X-< J-1 )*XSTEP+K*DELTA 

INDEX  *  ISHF-K  NPTS+1  >*J-NPTS+K 

KBUF< INDEX)  -  F< A,X, 

100  CONTINUE 

200  CONTINUE 

201  L  -  L+1 
CALL  DSPLAY 

IF< I0PT . EQ . 1 HN )  RETURN 
IF<  L . GT .479  )  RETURN 


0239  CALL  EXEC<  2 , LUD, KBUF, 512, I TK, I  SECT) 

0240  RETURN 

0241  END 

0242  C 

0243  c*************************  ******* 

0244  C 

0245  SUBROUTINE  SOLVE<I,J> 

0246  C 

0247  C  This  subroutine  solves  the  IDIM-2  simultaneous  equations 
0248  C  for  the  second  derivatives  determined  by  the  data  values 
0249  C  in  row  I  or  column'J  of  the  input  buffer  IBUF. 

0250  C 

0251  COMMON  /UNIG/NPTS , IDIM, CC0EF< 64, 4 >, DELTA, ISHF, XSTEP, 

0252  *  DERIV<64,62>,T<62,62> 

0253  COMMON  /INPUT/IBUF< 64, 64 > 

0254  COMMON  /RR0(J/A< 4  >, S< 62 >, Y< 62 > 

0255  C 

0256  LL  -  IDIM-2 

0257  IF< I . EQ . 0 )  GO  TO  200 

0258  C 

0259  C  Compute -second  derivatives  for  row  I  of  input  buffer. 

0260  C 

0261  DO  50  K-1,LL- 

0262  Y<K>  -  IBUF<  I , K+2 >-2*IBUF<  I,K+1  >+IBUF<  I,K> 

0263  Y<K>  -  Y<  K  >*6/<  XSTEP*XSTEP  > 

0264  50  CONTINUE 

0265  C 

0266  DO  100  K-1,LL 

0267  S<K>  -  0.0 

0260  DO  75  N-1,LL 

0269  S<K)  -  S< K  >*T<  K,  M  )*Y< M > 

0270  75  CONTINUE 

0271  100  CONTINUE 

0272  RETURN 

0273  C 

0274  C  Compute  second  derivatives  for  column  J  of  input  buffer. 

0275  C 

0276  ?00  CONTINUE 

0277  DO  250  K-1, LL 

0278  Y<  K  >  -  IBUF<K+2, J >-2*IBUF< K+ 1 , J>+IBUF<K, J> 

0279  Y<  K  >  ■  Y<K>*6AXSTEP*XSTEP> 

0280  250  CONTINUE 

0201  C 

0282  DO  300  K-1,LL 

0283  DERIV<  J,  K  )  -0.0 

0204  DO  275  M«1,LL 

0285  . DERI V<  J ,  K)  -  DERIV< J, K  )+T< K, M >*Y< M ) 

0286  275  CONTINUE 

0207  300  CONTINUE 

0288  RETURN 

0289  END 

0290  C 

0291  c******************************** 

0292  C 

0293  SUBROUTINE  GETC0<I,J,N> 

0294  C 

0295  C  This  subroutine  computes  coefficients  needed  to  interpolate 
0296  C  span  N  along  row  ,1  or  column  J. 

0297  C 

0298  COMMON  /UNI Q/NPTS, IDIM, CC0EF< 64 , 4 >, DELTA, ISHF, XSTEP, 
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0299  *  DERIV<64,62>,T<62,62> 

0300  COMMON  *'INPUT/IBUF<  64, 64  > 

0301  COMMON  /RR0W/A<4>,S<62>,Y<62> 

0302  C 

0303  IF< I . EQ . 0  >  GO  TO  100 

0304  C 

0305  C  Row  interpolation 
0306  C 

030?  A<  1  >  «  0.0 

0308  A<  2  >  -  0.0 

0309  IF<N.GT.1>  A<  1  >  -  S<N-1)/'6.0 

0310  IF<N.LT.< IDIM-1 >>  A<  2  >  »  S<NV6.0 

031 1  A<3)  -  IBUF<  I ,  N  > 

0312  A<4>  -  IBUF<  I ,  N+1  ) 

0313  RETURN 

0314  C 

0315  C  Column  interpolation 
0316  C 

0317  100  CONTINUE 

0318  CC0EF< J, 1  )  -  0.0 

0319  CC0EF<J,2>  -0.0 

0320  IF<  N  .  GT .  1  >  CCOEF<J,1>  -  DER I  V<  J,  N-1  V6 . 0 

0321  IF< N . LT .63 >  CC0EF<J,2>  -  DERI V< J, N >/6 . 0 

0322  CCOEF<  J, 3  >  -  IBUF<N,J> 

0323  CCOEF<  J, 4  >  -  IBUF<N+1,J> 

0324  RETURN 

0325  END 

0326  C 

0327  C******■*■*************♦***,4■*******,•,,•,* 

0328  C 

0329  FUNCTION  F<A,X, I  ) 

0330  C 

0331  C  Thi*  routine  .computes  the  interpolated  value  for  x  coordinate  X 
0332  C  which  lies  in  the  Ith  span. 

0333  C 

0334  COMMON  ✓UNIQ/NPTS, IDIM, CC0EF< 64, 4 >, DELTA, ISHF, X3TEP , 

0335  *  DERIV<  64, 62  >,  T<  62, 62  > 

0336  DIMENSION  A<4> 

0337  H  -  XSTEP 

0338  DX1  -  X-< 1-1 >*H 

0339  DX2  ■  H-DX1 

0340  F  *  A<  I  >*DX2*<DX2+DX2.'H-H> 

0341  *  +A<  2  >*DX 1 *<  DX 1 *DX 1 /H-H  > 

0342  +  +A<  3  >*DX2#/H+A<  4  )#DX  1  /H 

0343  IF<  F . LT . 0 . 0  >  F-0.0 

0344  RETURN 

0345  END 

0346  C 

0347  C**************************  ********* 

0348  C 

0349  SUBROUTINE  TREAD< ICNUD, I , IP  > 

0350  C 

0351  C  This  subroutine  reads  one  row  of  the  input  picture, 

0352  C  assuming  the  following  formati 

0353  C  1  row  -  32  16-bit  words  -  64  8-bit  pixels 

0354  C 

0355  COMMON  /'INPUT/IBUF<  64 , 64  > 

0356  COMMON  /UNIG/NPTS, IDIM, CC0EF<64, 4 >, DELTA, ISHF, XSTEP, 

0357  *  DERIV<64,62>,T<62,62> 

0358  DIMENSION  LBUF<  32  >, IR<  2  >, IP<  5  > 
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0359 

0360 

0361 

0362 

0363 

0364 

0365 

0366 

0367 

0363 

0369 

0370 

0371 

0372 

0373 

0374 

0375 

0376 

0377 

0378 

0379 

0380 

0381 

0382 

0383 

0384 

0385 

0386 

0387 

0388 

0389 

0390 

0391 

0392 

0393 

0394 

0395 

0396 

0397 

0398 

0399 

04  00 

0401 

04  02 

04  03 

04  04 

04  05 

0406 

04  07 

04  08 

04  09 

0410 

0411 

0412 

0413 

0414 

0415 

0416 

0417 

0418 


EQUIVALENCE  <IR,REC> 

REG  =  EXEC< 1 , ICNUD, LBUF, 32  > 

IF< IAND< IR, 2  00B  ) . EQ . 20  0B >  GO  TO  995 
LL  *  IDIM72 

IF< IR<  2  > . GE . LL  >  GO  .  TO  75 
N  »  IR<2>-H 
DO  50  K«N, LL 
50  LBUF< K >  »  0 
C 

75  CONTINUE 

DO  100  J-1,LL 

IBUF< I , 2* J-1 >  *  IAND< LBUF<  J >7256, 3?7B >*4 
IBUF< 1 , 2*  J  >  »  IAND<LBUF< J>,3?7B>*4 
100  CONTINUE 
RETURN 
C 

995  URITE< IP, 996 ) 

996  FORMAT<  "£t>a1 9rOC£KInval  id  #  of  records  in  picture  file") 

STOP 

END 

C 

c 

SUBROUTINE  DSPLAY 
C 

C  This  subroutine  nakes  the  appropriate  EXEC  call  to  write  a  line  of 
C  the  output  picture  to  the  selected  screen,  specified  by  NPRT  and  LU3. 
C 

COMMON  70UTPT7I0PT,LUD, IPIC,L, JBUF< 51 4 >, NPRT, LU3, ITK, ISECT 
DIMENSION  KBUF<  512) 

EQUIVALENCE  <  KBUF< 1 >, JBUF< 3 ) > 

C 

ITK  -  L71 2+< IPIC-1 >*40 
ISECT  -  <L-< ITK-< IPIC-1 >*40>*12>*8 
IF<  NPRT . EQ . 2HRA  >  GO  TO  200 
C 

CALL  EXEC<  2, LU3, KBUF, 51 2, L  > 

RETURN 

C 

200  CALL  EXEC<2,LU3, JBUF, 514) 

RETURN 

C^2 

END 

C 

C*********************************** 

c 

BLOCK  DATA 

COMMON  7INPUT7IBUF<  64, 64  > 

COMMON  70UTPT7I0PT, LUD, IPIC,L, 0BUF<  514), NPRT, LU3, ITK, ISECT 
COMMON  7UNIQ7NPTS, IDIM, CC0EF< 64,4  >, DELTA, ISHF,XSTEP, 

*  DERIV<  64, 62 ), T<  62, 62 ) 

COMMON  7RR0W7AC  4  ), S<  62 ), Y<  62  > 

OATA  NPTS737, IDIM7647,T70. 07,DELTA70.57 
DATA  NPRT72HRA7, LU371 2/ 

DATA  JBUF< 1  V5001 87, JBUF< 2  >7 1 024/ 


DATA 

T< 

1  , 

1 

>74 

.  07, 

T< 

1  , 

,2 

>71 , 07, 

* 

T< 

2, 

1 

>71 

.  07, 

T< 

2, 

,2 

>74. 07, 

T< 

2, 

3 

>71 

,  07 

* 

T< 

3, 

2 

>71 

.  07, 

T< 

3, 

,3 

>74, 07, 

T< 

3, 

4 

>71 

,  07 

* 

T< 

4, 

3 

>71 

.  07, 

T< 

4, 

,4 

>/A . 07, 

T< 

4, 

5 

>71 

,  0/ 

* 

T< 

5, 

4 

>71 

.  07, 

T< 

5, 

,5 

>74, 07, 

T< 

5, 

6 

>71  , 

,  07 
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T<6,5>71  .  07, T< 6, 6 >74, 07,  T<  6, 7  >71  .  07, 
T<7,6>71  .  07, T<  7, 7  >74 . 07,  T<  7,8)71  .  07, 

T<  8, 7)71  ,  07,  T<  8, 8  >74. 07,  T<  8, 9)71 . 07, 

T<  9, 1 >71 . 07, T< 9,9 >74. 07, T< 9, 10)71 .07, 

T<  1  0,9>71 . 07,  T<  1  0,  1  0>74. 07,  T<  1  0,  1  1  >71 . 07, 
T< 1 1 , 1 0>71 . 07, T< 1 1 , 1 1 >74. 07, T< 1 1 , 12)71 . 07, 
T< 12, 1 1 >71 . 07, T< 12, 12 >74. 07, T< 12, 13)71 . 07, 
T<  13, 12)71 ,07,T<  1 3, 13  >74.07,  T<  13,  14>71 .07, 
T<  14, 13)71 . 07,  T<  14,  14)74. 07,  T<  14, 13  >71 . 07, 
T<  15,  14)71 .07,  T<  15,  15)74. 07,  T<  15,  16)71 . 07, 
T<  16,  15)71 . 07,  T<  16,  16)74. 07,  T<  16,  17)71 . 07, 
T<  17,  16)71 . 07,T<  17,  17)74. 07,T<  17,  18)71 . 07, 
T<  18,  17)71 . 07,  T<  18,  13)74. 07,  T<  18,  19)71 . 07, 
T<  19,  18)71 .07,  T<  19,  19)74. 07,  T<  19,20)71 . 07, 
T<  20,  19)71 . 07,  T<  20, 20  >74. 07,  T<  20,  21  >71 . 07, 
T<21, 20)71 .07,T<2!,21  >74 . 07,  T<  21 , 22  >71 .07, 
T<  22,2 1 >71 . 07, T< 22, 22 >74. 07, T< 22, 23)71 . 07, 
T< 23, 22 >71 . 07, T< 23, 23 >74,07, T< 23, 24 >71 . 07, 
T< 24, 23)71 . 07, T< 24, 24 >74. 07, T< 24, 25 >71 . 07, 
T< 25, 24)71 . 07, T< 25, 25 >74. 07, T< 25, 26 >71 . 07, 
T< 26, 25)71 . 07, T< 26, 26 >74. 07, T< 26, 2? >71 . 07, 
T<  27,26  >71 . 07, T<27, 27)74, 07, T< 27, 28 >71 . 07, 
T< 28, 27)71 . 07, T< 28, 28 >74. 07, T< 28, 29 >71 . 0, , 
T< 29, 28)71 . 07, T< 29, 29 >74. 07, T< 29, 28 >71 . 07, 
T<  30, 29  >7 1  .  07,  TOO,  3  0)74. 07,  TOO,  31  >71 . 07, 
T<31 ,31  >71 . 07, T<  31 ,31  >74. 07,T<31 ,32  >71 . 07, 
T<  32, 3 1 >71 . 07, T<  32 , 32  >74 , 07, T<32, 33  >71 . 07, 
T< 33, 32  >7 1 . 07, T< 33, 33  >74 . 07, T<33, 34  >71 .07, 
T<  34,33  >71 . 07, T<  34, 34 >74 . 07, T<34, 35 >71 .07, 
T<  35, 34  >71 . 07,  T< 35, 35 >74. 07, T< 35, 36 >71 . 07, 
T<  36, 35  >71 . 07,  T06,  36  >74 . 07,  T06, 37  >71 . 07, 
T<  37 , 36  >71 . 07, TO?, 37  >74. 07, K37, 38  >71 .07, 
T08, 37  >71 . 07,  T<  38, 38  >74 . 07,  T<38, 39  >71 . 07, 
T<  39, 38  >71 . 07,  T<  39,  39  >74 .07,T<39,40>71 . 07, 
T< 40, 39 >71 . 07, T< 40, 40 >74. 07, T< 40, 41 >71 . 07, 
T<41 ,40)71 . 07,T<41 ,41  >74 . 07,  T<  41 , 42  >71 .07, 
T<  42, 4 1 >71 . 07, T< 42, 42)74. 07, T< 42,43)71 . 07, 
T< 43, 42 >71 . 07, T< 43, 43 >74. 07, T< 43, 44 >71 . 07, 
T< 44, 43 >71 . 07, T< 44, 44  >74 . 07, T< 44, 45 >71 . 07, 
T<45, 44)71 . 07, T< 45, 45 >74.07, T< 45, 46 >71 .07, 
T< 46, 45 >71 . 07, T<  46, 46 >74 . 07, T< 46, 47 >71 . 07, 
T< 47, 46 >71 . 07, T< 47, 47 >74 . 07, T<  47, 48  >71 . 07, 
T< 48 , 47  >71 . 07, T<  48,48  >74 . 07, T<  48, 49  >71 . 07, 
T<  49,48)71 . 07, T<49, 49)74. 07, T<  49,50)71 .07, 
T< 50, 49 >71 .07, T< 50, 50 >74. 07, T< 50, 51 >71 . 07, 
T<  51, 50  >71 .07,T<51,51  >74 . 07,  T<  51 , 52  >71 .07, 
T<  52, 5 1  >71 . 07,  T<  52, 52  >74. 07,  T<  52,  53  >71 . 07, 
T< 53, 52)71 . 07, T< 53, 53 >74. 07, T< 53, 54 >71 .07, 
T< 54, 53 >71 . 07, T< 54, 54 >74. 07, T< 54, 55 >71 . 07, 
T< 55, 54 >71 . 07, T< 55, 55 >74. 07, T< 55, 56 >71 . 07, 
T< 56, 55 >71 . 07, T< 56, 56 >74. 07, T< 56, 57 >71 . 07, 
T< 57,  56 >71 . 07,  T< 57, 57 >74. 07, T< 57, 58 >71 . 07, 
T< 58, 57 >71 . 07, T< 58, 58 >74. 07, T< 58, 59 >71 . 07, 
T< 59, 58 >71 . 07, T< 59, 59 >74. 07, T< 59, 60 >71 .07, 
T<  60, 59  >71 . 07,  T<  6  0,6  0  >74. 07,  T<  60, 61  >71 .07, 
T< 6 1,60 >71 ,07,T<61,61  >74 . 07, T< 61 , 62 >71 .07, 
T<  62, 6 1 >71 . 07, T< 62, 62 >74. 07 


-31- 


******  T  =  uGClOCl  IS  ON  LU  04 


0001 
0002 
0003 
0  0  04 
0005 
0006 
0007 
0006 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0  032 
0033 
D034 
D035 

3036 

3037 

3038 

3039 

3040 
0041 
0042 
0043 
0044 
0045 
0046 
C'047 
0048 
0049 
0050 
0051 
3052 
0053 
0054 

3055 

3056 

3057 

3058 
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c* 
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* 

c* 
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♦ 

c* 
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* 

c* 
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c*** ****************************************************** ********** 

C*  This  program  is  used  to  document  Program  EXPND  * 

C* ************************************************************ ****** 

c 

01  Software  Blueprint  —  Level  B  Design  of  the  Cubic  Spline  Expansion 
02  Program  Modules 

03  Module  Declaration 

#1  10  MAIN 
#2  20  INVRS 

#3  20  TREAD 

#4  20  SOLVE 

#5  20  CETCO 

#6  20  ROU 

#7  20  F 

#8  20  DSPLAY 

03  Module  Reference  Structure 

<calling  proc  name> f i -<cal led  proc  name  list> 

#1  10  MAIN i i "INVRS, DSPLAY, TREAD, SOLVE , ROW , CETCO, F 
42  20  INVRS 

#3  20  TREAD 

44  20  SOLVE 

45  20.GETCO 

46  20  ROW . i -SOLVE, CETCO, F, DSPLAY 

47  20  F 

48  20  DSPLAY 

02  Data . 

03  Data  Declaration  —  all  names  follow  FORTRAN  default  data  types 


main  program 

matrix  inversion 

read  one  row  of  input  tape 

solve  matrix  eqs  for  2nd  derivatives 

compute  interpolation  coefficients 

interpolate  across  row  of  data  points 

evaluate  interpolated  value 

output  interpolated  row  to  screen 


Buf fer 
DELTA, 
ICNWD, 
IDIM, 
INDEX, 
IPIC, 
ISECT, 
I  SHF, 

ITK, 

L, 

LDISK, 

LIM, 

LUD, 

LU3, 

MT, 

NPTS, 


relative  spacing  between  adjacent  output  pixels 
control  word  for  EXEC  calls  to  tape  drive  MT 
4pixels  per  side  of  input  picture 
pointer  into  output  array  KBUF 
output  picture  number 

sector  number  for  storing  picture  on  disc 

bias  count  needed  to  center  output  picture  row  in 

array  KBUF 

track  number  for  storing  picture  on  disc 
line  count  for  COMTAL  or  disc  output  control 
4  of  extra  rows  on  top  &  bottom  of  screen  when 
output  picture  is  centered 

4  of  extra  rows  on  top  &  bottom  of  disc  space 

when  output  picture  is  centered 

logical  unit  number  of  output  disc 

logical  unit  number  of  output  screen 

Initialize  LU3  to  12 

logical  unit  number  of  input  tape 

4  of  data  points  to  be  inserted  between  adjacent 


mm 

^  '  V  ■  ■ . .  ■  v 

\r rr v 

». '  r  T?  V  1  ■/M':*  V*  ■  T  T7  ».■  -j;  J  ■  i»f'  ,  ■ — »~r»  4,-  .■» 

1^' 

0059 

input  data  pixels 

WH- 

0060 

Initialize  NPTS  to  3 

0061 

X, 

position  coordinate  of  interpolated  pixel 

0062 

XSTEP, 

relative  spacing  between  adjacent  input  pixels 

0063 

0064 

Array 

0065 

A(4>, 

coefficients  needed  to  interpolate  between  2 

Urn 

0066 

pixels  in  a  row 

i 

006? 

CC0EF< IDIM, 4  >,  for  each  column,  coefficients  needed  to 

0068 

interpolate  between  2  pixels  in  that  column 

0069 

DERIV< IDIM, IDIM-2 ),  for  each  column,  2nd  derivatives  at  each 

0070 

input  pixel  (assuming  0  at  first  and  last 

0071 

pixels) 

0072 

IBUF< IOIM, IDIM >,  stores  entire  input  picture 

1  li 

0073 

IP<  5  >, 

stores  RMPAR  parameters  —  only  first  one  is 

(\ 

0074 

used;  IP  ■  logical  unit  #  of  user's  terminal 

• 

0075 

JBUF(514),  output  array  for  RAMTEK;  first  2  words  are 

•v 

0076 

predetermined,  rest  are  data  words 

0077 

Initialize  JBUF( 1 >  to  octal  5001  and  JBUF<2> 

0078 

to  decimal  1024 

g 

0079 

KBUF<512),  output  array  for  COMTAL  or  disc 

0080 

S< IDIM- 

*2),  2nd  derivatives  at  each  input  pixel  in  a  row 

» 

0081 

(assuming  0  at  first  and  last  pixels) 

0062 

T< IDIM- 

■2, IDIM-2)  inverse  of  matrix  arising  from  the  IDIM-2 

0083 

simultaneous  equations  for  the  2nd  derivativ 

0  084 

for  a  cubic  spline  fit  through  IDIM  equally 

0085 

spaced  data  points  with  2nd  derivative  of  0 

4 

0086 

at  the  endpoints 

^  ' 

008? 

Initialize  T  to  4  along  the  diagonal,  1 

0088 

off  the  diagonal,  and  0  elsewhere 

*  .' 

0069 

03  Data  Reference  Structure 

0090 

<data  structure  typeXdata  clement  name) i j *<ref erencing  proc's) 

0091 

m 

0  092 

Buffer 

0093 

DELTA  j 

i*  MAIN, ROU 

0094 

ICNUD  : 

!«  MAIN, TREAD 

0  095 

IDIM  i 

i«  MAIN, TREAD, SOLVE, GETCO, ROU 

'• 

0096 

INDEX  s 

i-  MAIN, ROU 

009? 

IPIC  i 

s=  MAIN, DSPLAY 

a 

0098 

ISECT  i 

:  =  MAIN, ROU, DSPLAV 

m 

0099 

I  SHF  : 

»-  MAIN, ROW 

y. 

01  00 

ITK  ■ 

i«  MAIN, ROU, DSPLAY 

*„• 

01  01 

L  i 

i-  MAIN, ROU, DSPLAY 

01  02 

LDISK  i 

:»  MAIN 

Ot  03 

LIM  i 

i-  MAIN 

01  04 

LUD  i 

i ■ ■ MA IN, ROU 

j 

01  05 

LU3  i 

s-  MAIN, DSPLAY 

01  06 

MT  i 

«»  MAIN 

01  07 

NPTS  i 

i-  MAIN, ROU 

01  08 

X  : 

i-  MAIN, ROU, F 

01  09 

XSTEP  i 

i-  MAIN, SOLVE, ROW, F 

01  1  0 

i 

0111 

Array 

■.-; 

0112 

A  : 

:■  MAIN, GETCO, ROW, F 

• 

01  13 

CCOEF  » 

i-  MAIN, GETCO 

■.- 

01  14 

DERIV  t 

i»  SOLVE, GETCO 

*  r 

0115 

I8UF  s 

i-  MAIN, TREAD, SOLVE, GETCO 

■  •; 

01  (6 

IP  i 

i-  MAIN, CHECK, TREAD 

4 

01  1? 

JBUF  i 

i-  DSPLAY 

01  18 

KBUF  i 

i-  MAIN, ROW 

,•*■ 
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01  19 
0120 
0121 
0122 
0123 
0124 
0125 
0126 
012? 
0126 
0129 
0130 
0131 
0132 
0133 
0134 
0135 
0136 
0137 
0139 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 
0161 
0162 
0163 
0164 
0165 
0166 
0167 
0168 
0169 
0170 
0171 
0172 
0173 
0174 
0175 
0176 
0177 
0178 


S  m»  SOLVE, GETCO 

T  ii-  INVRS, SOLVE 

02  Control 

03  Control  Declaration 
Switch 

IOPT  Of  Status  < "V" ,  "N" >  "save  on  disc"  response 
NPRT  Of  Status  <"RA" , "CO" >  output  screen  identifier 
03  Contro1  Reference  Structure 

<switci.  name>  i  i  «<where  -  set >/< where  tested'/ 

IOPT  li-  MAIN/MAIN, ROW 
NPRT  i  i-  MAINXMAIN, DSPLAY 
02  Procedure . Definition 
#1  Procedure  MAIN 

Questions  asked  .of  user  at  run  timet 

<1>  "program  .defaults  to  RANTEK  screen  on  LU  12. 

Do  you  want  to  change  either  default?  Y/N" 

If  answer  is  "V",  <1a>  "input  CO  for  COMTAL  or  RA  for  RAMTEK" 

< 1 b  >  "input  LU  #" 

<2>  "specify  expansion  factor  —  e.g.  2  will  approximately  doubl 
input  picture  which  is  64x64.  < or  just  RETURN  and  use  defat 
of  4" 

<  3  >  "specify  input  picture  tape  unit" 

<4>  “do  you  want  to  save  output  picture?  Y/N" 

If  answer  is  "Y",  <4a>  "specify  output  disk  LU  and  picture  numb« 
At  completion  of  expansion, 

<5>  "do  you  want  to  expand  the  next  picture?  Y/N" 

Call  RMPARj  get  LU  #  of  user's  terminal 

Write  questions;  prompt  user  for  specifics 

Read  answers; 

Call  INVRS;  invert  matrix  T 

Initialize  L, ISHF,XSTEP>  and  ICNWD; 

If  NPRT  »  "RA" 

Then  issue  reset  command  to  RAMTEK; 

End  If; 

Write  leading  rows  of  0  to  screen; 

If  IOPT  »  “Y" 

Then  write  leading  rows  of  0  to  disc; 

End  If; 

Loop  until  entire  input  picture  has  been  read 

Call  TREAD;  read  next  row  of  input  picture 

End  Loop; 

Loop  until  all  columns  of  input  picture  have  been  processed 
Call  SOLVE;  get  2nd  derivatives  at  column  pixels 

< 0  at  endpoints) 

End  Loop; 

Embed  1st  row  of  input  picture  in  output  array  K8UF, 
leaving  NPTS  gaps  in  between  for  interpolated  values; 

Call  ROW;  fill  inthis  row 

Loop  until  all  input  rows  have  been  processed 

Loop  until  all  input  columns  have  been  processed 

Call  GETCO;  get  coefficients  to  interpolate  along  this 

column  between  old  and  current  input  rows 

End  Loop; 

If  NPTS  is  not  0  Then 

Loop  until  NPTS  pixel  positions  have  been  processed 
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X  i-  X+DELTA;  put  X  at  next  pixel  position  j 

Loop  until  all  input  columns  have  been  processed 

Call  F;  interpolate  at  X  in  this  column  | 

Urite  new  value  in  output  array  KBUF; 

Replace  the  data  point  in  the  input  array  1BUF  at 
this  column  and  old  row  by  the  new  value; 

End  Loop; 

Call  ROW;  interpolate  across  the  row  jus‘1  manuf  actured; 

Move  current  row  pointer  into  old  row  pointer; 

End  Loop; 

End  If; 

Embed  current  input  row  in  output  array  KBUF, leaving 
NPTS  gaps  in  between  for  interpolated  values; 

Call  ROW;  fill  in  this  row 

End  Loop; 

Urite  trailing  rows  of  0  to  screen; 

If  IOPT  -  "V" 

Then  write  trailing  rows  of  0  to  disc; 

End  If; 

End  MAIN;  , 

#2  Procedure  INVRS 

This  subroutine  was  obtained  verbatim  from  NUMERICAL  METHODS 
by  Robert  U.  Hornbeck . 

This  subroutine  is  only  called  at  initialization  to  invert 
matrix  T. 

Array 

A< IDIM-2, IDIM-2  ),  original  matrix 
J< IDIM-2+21 >  temporary  storage 

Buffer 

M  dimension  of  matrix  to  be  inverted 

Calculate  .inverse  .of  A; 

Store  inverse  in  A; 

Return; 

End  INVRS; 

#3  Procedure  TREAD 

Array 

LBUF< I  DIM/2 )  stores  one  row  of  input  picture 

Read  next  record  from  tape  MT ; 

If  EOF  encountered 

Then  write  "invalid  #  of  records  in  picture  file"; 

Stop; 

End  If; 

If  #  of  words  read  <  32 

Then  zero-fill  rest  of  input  array; 

End  If; 

Repack  each  byte  in  LBUF  into  a  word  in  IBUF,  preserving 
the  order  of  the  bits; 

Shift  each  word  left  2  bits  to  rescale  data  to  10  bits; 

End  TREAD; 

#4  Procedure  SOLVE 


Let  g<x>  be  the  cubic  spline  fit  through  an  equally  spaced 
line  of  pixels  located  at  x  ,x  ,...,x  and  with  gray  levels 


0239  G  ,G  ,  ...,G  .  Then  the  2nd  derivatives  of  g  at  the  pixels 

0240  1  2  N 

0241  satisfy  the  following  equations: 

0242  g"<x  >+4+g"<x  >+gM<x  )  *  6*E<G  -2*G  +G  >A'x>+*23 

0243  k  k  +  1  k+2  k+2  k  +  1  k 

0244 

0245  for  k*1 , 2, , . . , N-2 .  Note:  'x  »  delta  x. 

0246  There  arc  N-2  equations  in  the  H  unknowns  g"<x  >,k*1 , 2,  , .  , ,  N . 

0247  k 

0248  The  additional  constraints  of  2nd  derivative  •  0  at  the  end- 

0249  points  permit  the  solution  of  these  equations.  In  matrix  form, 

0250  the  equations  are: 

0251  <4  I  0  .  .  .  OX  g"<2>  > 

0252  <\  4  f  0  .  OX  g"<3>  ) 

0253  <0  1  4  1  0..X  g"<4>  >  =  Y,  Y<  k  >  »  6+C<G  -2*G  +G  >/<Ax>++23 

0254  < . IX . >  k+2  k  +  1  k 

0255  <0  ...  0  1  4  X  g“<  N- 1  >> 

0256 

0257  In  this  application,  N  *  IDIM  and  T  is  the  inverse  of  the 

0258  coefficient  matrix. 

0259  Buffer 

0260  I,  0  or  #  of  the  row  in  input  array  IBUF 

0261  containing  the. desired  pixels 

0262  J  0  or  #  of  the  column  in  input  array  IBUF 

0263  containing  the  desired  pixels 

0264  Array 

0265  Y<IDIM-2>  temporary  storage  for  expressions  on  right 

0266  hand  side  of  equations  for  g 

0267  If  I  is  non-zero  row  interpolation 

0268  Then  compute  Y ;  note  G  =  IBUF<I,K>  and  Ax  =*  XSTEP 

0269  S  T+Y;  k 

0270  Return; 

0271  Else  compute  Y;  column  interpolation 

0272  note  G  -  IBUFCK,J>  and  Ax  =  XSTEP 

0273  k 

0274  DERIV< J, . )  =  T*Y; 

0275  Return; 

0276  End  If; 

0277  End  SOLVE; 

0278  #5  Procedure  GETCO 

0279  g"  has  been  determined  at  the  pixel  values  along  either  a  row 

0280  or  column  of  the  p i ctur e  and  stored.  This  de term i nes  g<  x  >  on 

0281  each  span  (between  consecutive  input  pixels)  because  each  cubic 

0282  can  be  written  as  a  function. of  g": 

0283  between  pixels  :at  x  and  x  < span  k >, 

0284 

0185  k  k+1 

0286  Note:  Ax  *  delta  x 

0287  g<x)  ■  g"<x  )/6*C<x  -x>**3/Ax  -  Ax+<x  -x>3  + 

0288  k  k+1  k+1 

0289 

0290  g"<  x  >/6*t<  x-x  >*+3  7Ax  -  Ax*<x-x  >3  + 

0291  k+1  k  k 

0292 

0293  G  +C< x  -xVAx3  +  G  *C<x-x  >/Ax3 

0294  k  k+1  k+1  k 

0295  Buffer 

0296  I,  either  0  or  #  of  row  in  input  array  IBUF 

0297  containing  the  desired  pixels 

0298  <i>  either  0  or  #  of  column  in  input  array  IBUF 


0299 

0300 

0301 

03  02 

0303 

0304 

0305 

0306 

030? 

0308 

0309 

031  0 

031  J 

0312 

0313 

0314 

0315 

0316 

0317 

0318 

0319 

0320 

0321 

0322 

0323 

0324 

0325 

0326 

0327 

0328 

0329 

0330 

0331 

0332 

0333 

0334 

0335 

0336 

0337 

0338 

0339 

Q340 

0341 

0342 

0343 

0344 

0345 

0346 

074? 

0348 

0349 

0350 

0351 

0352 

0353 

0354 

0355 

0356 

0357 

0358 


containing  the  desired  pixels 
N  span  11 

If  I . is  non-zero 

Then  row  interpolation 

If  N  >  1  note  g"  is  0  at  x  and  x 

1  IDIM 

Then  A<  1  >  :»  S<H-1  >/6; 

Else  A<  1  >  i *  0) 

End  If) 

If  N  <  IDIM-1 

Then  A<  2 )  i-  S<N>/6; 

Else  A<2>  :a  0; 

End  If) 

A<  3  >  i-  IBUF<  I,N>j 
A<  4  )  :=  IBUF<  I  >  N+1  >; 

Return; 

Else  column  interpolation 

If  N  >  1  note  g"  is  0  at  x  and  x 

1  IDIM 

Then  CC0EF< J, 1 >  !-  DERI V< J, N-1 V6; 

Else  CCOEF<  J,  1  )  i  -  0; 

End  If; 

If  N  <  IDIM-1 

Then  CC0EF<J,2>  :=  DERI  V<  J ,  N  )/6j 
Else  CCOEF< J,2  >  ! =  0; 

End  If) 

CC0EF<J,3>  j*  IBUF<N,J>; 

CC0EF<  J  j  4  >  !-  IBUF<N+1,J>) 

Return) 

End  If; 

End  GETCOj 
#6  Procedure  ROW 

NPTS  points  must  be  ,  inserted  between  each  pair  of  pixels  in 
the  designated  row  of  input  array  IBUF.  The  position  values 
assigned  to  these  pixels  arei 

Kth  pixel  has  x  coordinate  <K-1>*XSTEP. 

This  allows  NPTS  pixels  to  be  inserted  at  subintervals  DELTA 
in  each  XSTEP  interval. 

Buffer 

I  #  of  row  in  input  array  IBUF  containing 

the  desired  pixels 

If  NPTS  non-zero  bypass  processing  if  0  expansion  is  opted 

Then 

Call  SOLVE;  get  2nd  derivatives  for  cubic  spline  fit 

across  this  row 

Loop  until  last  pair  of  pixels  processed 

Call  CETCO)  get  the  4  coefficients  which  specify  the 

cubic  through  current  pair  of  pixels 
Loop  until  NPTS  points  have  been  inserted 
X  j-  X+DELTA; 

Call  F;  evaluate  cubic  at  X 

Store  in  output  . array  KBUF; 

End  Loop; 

End  Loop; 

End  If; 

Call  DSPLAY)  output  new  row  to  screen 

If  I0PT  -  "N" 
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Then  Return; 

Else  Store  row  on  disc  LUO; 

Return; 

End  If; 

End  ROLI; 

#7  Procedure  F 

evaluate  the  cubic  interpolation  polynomial 
Buffer 

X,  position  of  new  pixels  to  be  interpolated 

I ,  span  # 

0X1,  temporary  storage  of  distance  between  new 

pixel  and  left  endpoint  of  span 

0X2,  temporary  storage  of  distance  between  right 

endpoint  of  span  and  new  pixel 
H  temporary  storage  of  XSTEP 

Array 

A<  4  >  weighting  coefficients  for  this  span 

DX1  J*  X-< 1-1 >*H; 

0X2  i»  H-DX1 ; 

F  :*  A<  f  >*DX2*CDX2*+2/'H  -  HI  + 

A<2 >*DX1*CDX1**27H  -  HI  + 

A<  3  >*DX2/H  +  A<4>+DX1/H; 

Return; 

End  F; 

#8  Procedure  OSPLAY 

Compute  track  and  sector  for  next  row; 

If  NPRT  «  “RA” 

Then  output  JBUF  to  RAMTEK; 

Else  output  K8UF  to  CGMTAL; 

End  If; 

Return; 

End  OSPLAY; 
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0001 
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0003 

c* 

Uritten  by  Ginny  Kalb 

* 

0004 

c+ 

Technology  U.  S.  A. 

* 

0005 

c+ 

P.  0.  Box  55333  Ph  <301)  292-2592 

* 

0006 

c* 

Oxon  Hill  Station 

'¥ 

0007 

c+ 

Washington  D.  C.  20022 

* 

0008 

c* 

Rev .  29  Jan  1981 

0009  c*'***********************************1*****************************-*'* 

0010  C+  This  program  is  used  to  explain  how  to  use  Program  EXPND.  * 

00(1  c*** ********************************************** ******** ********** 

0012  C 

0013  User's. Guide  to  Cubic  Spline  Expansion 

0014  Program  EXPND 

00(5  This. program  expands. a  64x64  pixel  image  by  row  and  column  cubic  spline 
0016  interpolation.  The  expansion  factor  is  selectable,  ranging  from  1  (no  ex- 
0017  pansion)  to  8  <output  image  is  505x505  pixels).  The  output  picture  is 
0018  displayed  on  cither  a  RAMTEK  or  C0MTAL  screen  and  optionally  stored  on  dig 
0019  The  user  options  are  specified  at  run  time  in  response  to  prompts  issued 
0020  by  the  . program.  The  input  must  be  on  tap*  and  the  tape  must  be  positioned 
0021  by  the  user  to  the  desired  picture  file  prior  to  program  execution. 

0022  It  is  possible  to  speed  up  program  execution  if  the  user  has  a  data  tap( 
0023  with  consecutive  images  all  of  which  are  to  be  expanded  by  the  same  facto) 
0024  because  then  the  initial ization  of  the  algorithm  can  be  skipped. 

0025  Note)  the  input  tape  is  not  rewound  at  the  end  of  the  program. 

0026  Sample  computer-user  dialogue) 

0027  C  NOTE  program  defaults  to  RAMTEK  screen  on  LU  12. 

0028  Do  you  want  to  change  either  default?  Y/N 

0029 

0030  U  Y 
0031 

0032  C  .  input  CO  for  CONTAL  or  RA  for  RAMTEK 
0033 

0034  U  RA 
0035 

0036  C  . input  LU  # 

0  037 

0038  U  16 
0039 

0040  C  .specify  expansion  factor  —  e.g.  2  will  double  input  picture  which 
0041  is  64x64  < or  just  RETURN  and  use  the  default  of  4  ) 

0042 

0043  U  2 
0044 

0045  C  .specify  input  picture  tape  unit 
0046 

0047  U  7 

0048  delay  while  initialization  is  performed 
0049  C  .do  you  want  to  .save  output  picture?  Y/N 
0050 

0051  U  Y 
0052 

0053  C  .specify  output  disc  LU  and  picture  number 
0054 

0055  U  13,2 

0056  delay  while  expansion  is  performed 
0057  C  it  is. finished 

0058  .do  you  want  to  expand  the  next  picture?  Y/N 
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